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Abstract. We revisit the topic of the mean field ground state of a spin-1 atomic 
condensate inside a uniform magnetic field (B) under the constraints that both the total 
number of atoms (N) and the magnetization {A4) are conserved. In the presence of an 
internal state (spin component) independent trap, we also investigate the dependence 
of the so-called single spatial mode approximation (SMA) on the magnitude of the 
magnetic field and M. Our result indicate that the quadratic Zeeman effect is an 
important factor in balancing the mean field energy from elastic atom-atom collisions 
that are known to conserve both N and A4 . 



Mean field ground state of a spin-1 condensate in a magnetic field 



2 



1. Introduction 



Atomic Bose-Einstein condensates (BEC) liave provided a successful testing ground for 
theoretical studies of quantum many-body systems Ij. In earlier BEC experiments, 
atoms were spatially confined with magnetic traps, which essentially freeze the atomic 
internal degrees of freedom Most studies were thus focused on scalar models, 
i.e. single component quantum degenerate gases jS]. More recently, the emergence 
of spin-1 condensates [H El E] (of atoms with hyperfine quantum number F = 1) 
has created opportunities for understanding degenerate gases with internal degrees of 
freedom [3 |H1 El CHI HI El • 

In this paper, we investigate the mean field ground state structures of a spin-1 
atomic condensate in the presence of an external magnetic field (B). We focus on 
several aspects of the ground state properties strongly affected by the requirement 
that elastic atom-atom collisions conserve both the total number of atoms (A^) and 
the magnetization (Ai). Several earlier studies have focused on the global ground 
state structures when the conservation of Ai was ignored, or in the limiting case of 
a vanishingly small magnetic field (5 = 0) (HI El El HHl HH Hll- As we show in 
this study, in the presence of a nonzero magnetic field, the conservation of M. leads 
to ground state population distributions significantly different from that of the global 
ground state. 

Our system is described by the Hamiltonian (repeated indices are summed) [71 
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where ipjif) is the field operator that annihilates an atom in the j-th (j = +,0, 



internal state at location r, C: 



-rvV2M + V, 



cxt l 



6ij with M the mass of each 



atom and V^xt(^) an internal state independent trap potential. Terms with coefficients 
Co and C2 of Eq. (jH) describe elastic collisions of the spin-1 atom (|F = l,Mp = 
-|-,0, — )), expressed in terms of the scattering length Oq (02) for two spin-1 atoms in 
the combined symmetric channel of total spin (2), Co = 4:7rh^{ao + 2a2)/3M and 
C2 = 47r/i^(a2 — ao)/3M. Fri=x,y,z are spin-1 matrices with 
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The external magnetic field B is taken to be along the quantization axis 
induces a Zeeman shift on each atom given by 

f E+ \ 
Eo 


According to the Breit-Rabi formula ^21, the individual level shift can be expressed as 
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gi^iiB - -EhfsVI + a + a^. 
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i^o = - ^ - IerfsViT^, 

o Z 

E^= - ^ + gj^iiB - \EnYs^l-a + a\ (2) 

o Z 

where -Erfs is the hyperfine sphtting ^j, and gi is the Lande (yf-factor for the atomic 
nuclei with nuclear spin I. fij is the nuclear magneton and a = {gifijB + gjfisB) / E-ups 
with gj the Lande (^-factor for the valence electron with total angular momentum J. fiB 
is the Bohr magneton. 



2. Mean field approximation 

At near zero temperatures and when the total number of condensed atoms is large, the 
ground state is essentially determined by the mean field term $j = (ipi). Neglecting all 
quantum fluctuations we arrive at the mean field energy functional from Eq. [3 

Him] =Hs + EoN + - r/o(F,) + S{F^), (3) 



where the symmetric part 
H. = fdr 



(4) 



is invariant under the exchange of spin component indices, thus is independent of the 
external B field. The Zeeman shift as given by the Breit-Rabi formula ^ can be 
described by two positive parameters ^3] 

2r]o = E^- E+, 

25 =E+ + E.- 2Eq, (5) 

which measure approximately the linear and quadratic Zeeman effects. The -B-field 
dependence of tjq and S for a ®^Rb atom are displayed in Fig. [T] 

The elastic atomic collisions as described by the cq and C2 parts of the Hamiltonian 
(0) conserve both and Ai, which in the mean field approximation are given by 

N = J2 fdf{^p]{f)i,^{r))^ Yl ldr\^,{f^\\ 
M = I rff [(^i(f)^+(f)) - (^/'i(r)^/^_(f))] 



j dr[\^^{r)\'-\^^{r)\\ (6) 



Before continuing our discussion of the ground state structures, we shall first briefly 
comment on the importance of the above two constraints. In a typical experiment, the 
last stage before condensation consists of atomic evaporations, during which neither 
nor M. is conserved. For a scalar condensate, typically the ground state is obtained from 
a minimization of Eq. (jS)) subjected to the constraint of only A^ conservation. This gives 
rise to the Gross-Pitaevskii equation (GPE) and the associated condensate chemical 
potential, which mathematically is simply the Lagrange multiplier of the constrained 
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Figure 1. Approximate linear and quadratic Zeeman effects as characterized by 
parameters rjQ and S versus magnetic field B for a *^Rb atom. 

minimization. A spin-1 condensate requires the introduction of two Lagrange multipliers 
during the minimization subjected to both the and Ai conservation constraints, as 
was first performed in fH^ . 

When atomic interactions are ferromagnetic (c2 < as for ^''Rb atoms) and when 
the external -B-field is negligible, we have shown previously that the ground state 
structure is simply a state where all individual atomic spins are aligned in the same 
direction ^3]. In this case, the conservation of Ai can be simply satisfied by tilting 
the quantization axis away from the direction of the condensate spin. This can always 
be done if a system described by is rotationally symmetric, and thus contains the 
SO (3) symmetry [Zj. The presence of a nonzero B field, on the other hand, breaks 
the rotational symmetry, [e.g. the linear Zeeman shift, reduces the S0(3) to S0(2) 
symmetry] , thus the conservation of Ai has to be included in the minimization process 
directly. 

The global ground state phase diagram including both linear and quadratic Zeeman 
effect was first investigated by Stenger et al p^. In this early study, although the 
At conservation was included in their formulation, it was not separately discussed, 
consequently their results do not easily apply to systems with fixed values of Ai. The 
ground state structures as given in Ref. [HI correspond to the actual ground state as 
realized through a Ai non-conserving evaporation process (e.g. in the presence of a 
nonzero i?-field) that serves as a reservoir for condensate magnetization. Our study to 
be presented here, on the other hand, would explicitly discuss the phase diagram for 
fixed values of A4, which could physically correspond to experimental ground states 
(with/without a -B-field) due to a conserving evaporation process. Although more 
limited, as our results can be traced to linear trajectories of Ai = const, in the phase 
diagram of Ref. we expect them to be useful, especially in predicting ground state 
structures when a ready-made spinor condensate is subjected to external manipulations 



Mean field ground state of a spin-1 condensate in a magnetic field 5 

that conserve both and Ai. 

When atomic interactions are anti-ferromagnetic (c2 > 0), the global ground state 
was first determined to be a total spin singlet [Hj. More elaborate studies, including 
quantum fluctuations, were performed by Ho and Yip [T7] and Koashi and Ueda [THj . 
Unfortunately, these results jT3|THI do not correspond to actual ground states as realized 
in current experiments, because of the presence of background magnetic fields. For 
instance, the states as found in Ref. ^Tj are only possible if the magnetic field B 
is less than 70/iG at the condensate density as realized in the MIT experiments 
The Zeeman shift (see Fig. ^ due to the presence of even a small magnetic field can 
overwhelm atomic mean field interaction and typical atomic thermal energy, thus if it 
were not for the conservation of Ai, the ground state would simply correspond to all 
atoms condense into the lowest Zeeman sublevel of \Mp = 1). 

We now minimize H of Eq. ^ by denoting ^j{r) = ^jNj(f)j{r)e^^^\ with a real 
mode function 4>j{r) (/ (f)'j{r)df = 1) and phase 9j. It is easy to check that the phase 
convention of ferromagnetic/anti- ferromagnetic interactions as obtained previously ^Hj 
in the absence of a 5-field still remains true, i.e. 

9+ + 9- - 200 = 0, C2 < (ferromagnetic), (7) 
9+ + 9_ — 29q = vr, C2 > (anti- ferromagnetic). (8) 



3. Ground state in a homogeneous system 



In a homogeneous system such as a box type trap (of volume V), adopting the above 
phase convention, the resulting ground state energy functional becomes (+/— for C2 < 
and C2 > respectively) 

C2 



H[{N,}] = Hs + EoN + 



2V 



(iV+ - iV_)2 + 2No{^N+ ± ^JN^y 

-r7o(iV+-iV_) + 5(iV+ + iV_). (9) 

Expressing everything in terms of fractional populations and fractional magnetization 
Hi = Ni/N and m = M./N, and note that n+ + n_ = 1 — no, n+ — n_ = m. Equation 
(0) becomes 
H[{n,}] Hs 



N 



N 



+ Eq + c (n_|_ — + 2nQ{y/n^ ± y/nZY 

- ?7o(n+ - n_) + (5(n+ + n_), (10) 

with an interaction coefficient c = C2N/2V, tunable through a change of condensate 
density. 

We now minimize Eq. fllOj) under the two constraints n+ + no + = 1 and 
n+ — n_ = m. We restrict our discussion to the region — 1 < m < 1 as the special cases 
of m = ±1 are trivial. Because Hs, Eq, c, tjq, and m are all constants for given values 
of B, N, and V, the only part left to be minimized is 
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In the special case of c = 0, Equation (fTTj) reduces to 

T = 5{n+ + n_). (12) 

The ground state is then very simple. When 6 > 0, which seems to be always the case 
for quadratic Zeeman shift, the minimum is reached by having as large a Uq (thus as 
small a + n_) as possible, namely 

I I f I'^L rn>0 fO, m>0 

un = 1 — m , 77,+ = < , n_ = < , , j13 

' ' ^ \ 0, m < \ |m|, m < ^ ^ 

When 5 = 0, we have (in general) three condensate components with n± = (1— no±m)/2 
and < no < 1 — |m|. 

For ferromagnetic interactions with c < 0, we define x = + n_. The ground 
state is then determined by the minimum of 

J^ = g+{x)+6x, (14) 



with g+{x) = 2c(l — x){x + \/ x^ ~ m^). When 5 = 0, we find 

n± = ^(l±m)^ no = i(l-m2), (15) 

which is the same as obtained in |1H ITij . However with a nonzero 5 > 0, we find in 
general 

n± = ^(xo±m), no = 1 -xo > ^(1 - m^), (16) 

with xq being the root of equation g%[x) + 5 = 0, it turns out that there always exists 
one and only one solution to the equation. The equilibrium value for no is larger than 
the result of Eq. (fT^ because the quadratic Zeeman effect causes a lowering of the 
total energy if two \Mp = 0) atoms are created when an \Mf = +1) atom collides 
with an \Mp = —1) atom. Figure |21 displays the results of Eq. (fT^ for a typical 
®^Rb condensate, for which the atomic parameters are -Ehfs = (27r)6. 8347GHz jl2j . 

= 101. Sa^, and 02 = lOOAas [as is the Bohr radius) [15J. At weak magnetic fields, 
typically a condensate contains all three spin components. With the increasing of in- 
field, the quadratic Zeeman effect becomes important which energetically favors the 
|0) component, so typically only two components survive: the |0) component and the 
larger (initial population) of the |+) or |— ) component, so the ground state becomes 
(for m > 0) n+ ^ m and no — 1 — m. 

Finally we consider the case of anti-ferromagnetic interactions for c > 0, we have 

then 

J^ = g_{x)+6x, (17) 



with g-{x) = 2c(l — x){x — ya? — n?). For 5 = 0, we again recover the standard result 
no = 0, n± = ^(l±m), (18) 

if m 7^ 0. When m = 0, the ground state is under-determined as many solutions are 
allowed as along as they satisfy n+ = n_ = 1 — no with no G [0, 1]. 
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Figure 2. The dependence of fractional population for different spin component on 
TO and B for a spin-1 ^^Rb homogeneous condensate with N/V = 5 x lO^^cm"'^. 




Figure 3. The same as in Fig. 2, but now for a spin-1 Na condensate. 



In an external 5-field when 5 > 0, we first consider the special case of m = 0. 
It can be easily seen from Eq. (fT7|) that = 1 is the ground state. For m 7^ 0, we 
obtain the following result: when 6 > 2c[l — y/l^^mP], the ground state will have three 
condensate components with 

n± = (xo±m)/2, uq = 1 - xq, (19) 

where Xq is the root of equation g'_{x) -|- 5 = 0; When 5 < 2c[l — ^/l — m?], only |-|-) 
and |— ) components exist, i.e., n± = {1 ± m)/2. 

Figure El is the typical results for a spin-1 ^'^Na condensate. The atomic parameters 
are -Erfs = (27r)l. 7716GHz [T2], ao = SOa^ and 02 = SSa^ [I^. At 5 = there are 
only two condensate components, |-|-) and |— ). For B > Q but not very strong, there 
are two possibilities: three nonzero condensate components if m < and two nonzero 
condensate components if m > m-c, with 5{B) = 2c(l — y^l — ) . When 5-field gets 
stronger, i.e. S{B) > 2c, there are always three condensate components. 

Figure |3] summarizes the ground state structures of a homogeneous spin-1 
condensate in a -B-field for different c and m. 
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Figure 4. The ground state phase diagram for a homogeneous spin-1 condensate. 
Dashed curves and hues denote gradual transitions across the boundaries, sohd hues 
denote discontinuous jumps, xq is the sohition to equation g'^{x)~\-5 = and the curves 
for c > is determined by 5{B) = 2c[l — ^/l — m?]. The open circle at i? = 0, m = 
for c > denotes the family of degenerate ground state (-!-j^,no, ^~2° )- 



4. Ground state inside a harmonic trap 

In the previous section, we investigated in detail mean field ground state structures for a 
spin-1 condensate in a homogeneous confinement. For the case of a harmonic trap as in 
most experiments, there is no reason to believe a priori that the above conclusions still 
hold. In fact, the structures and phase diagrams as discussed before is only meaningful 
if the spatial mode function 0j(r) for different spin components is identical. Otherwise, 
it would be impossible to classify the rich variety of possible solutions. When the spatial 
mode functions are the same, the spatial confinement simply introduces an average over 
the inhomogeneous density profile of the mode function. 

The aim of this section, is therefore to determine the validity of the single mode 
approximation (SMA) in the presence of an external i?-field and a harmonic trap. For 
simplicity, we assume the trap to be spherically symmetric. We employ numerical 
methods to directly find the ground state solutions from the coupled Gross-Pitaevskii 
equation 

d 
d 

i^— $0 =[n + EQ + C2(n+ + n_)]$o + 2c2$o'^+$-, 
d 

z^— $_ = [n + E_+r] + C2{n_ + - + 02$^$% (20) 

subjected to the conservations of both and M [Eqs. Ti = -Ti^V"^ /2M + Vt{f) + 
Cqu, Uj = Vt{r) = Mu;^r^/2, and n = n+ + uq + n^. 77 is the Lagrange multiplier 

introduced to numerically enable the conservation of A4. 

It was shown previously that in the absence of an external S-field, and for 
ferromagnetic interactions, the SMA is rigorously valid despite the presence of a 
harmonic trap ^3]. We can also show that in the presence of a nonzero i?-field, the 
linear Zeeman shift does not affect the validity of the SMA because it can be simply 
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Figure 5. Typical densities of spatial mode functions for each components of a Wo 
(a) and a ^'^Na (b) condensate. The solid line denotes the |+) component, the dashed 
line the |— ) component, and the dash-dotted line the |0) component. The parameters 
are, N = 10^,^ = 2tt x 100 (Hz), B = 1.0 (Gauss), and m = 0.5. 

balanced by the external Lagrange multiplier t]. The quadratic Zeeman effect, on the 
other hand, can not be simply balanced, as it favors the production of two |0) atoms 
by annihilating one 1+) and one |— ) atom during a collision. Such unbalanced elastic 
collisions thus break the SO (3) symmetry of the freedom for an arbitrary quantization 
axis. Therefore, we do not in general expect the SMA to remain valid inside a nonzero 
5-field. 

Numerically, we find the ground state solutions of Eq. (pUj) by propagating the 
equations in imaginary time. We typically start with an initial wave function as that of 
a complex Gaussian with a constant velocity: exp[— (x^/2g^ + |/^/2g^ + z^/2g^) —ik-f]. 
Qx, %, Qz, and k are adjustable parameters which are checked to ensure that their choices 
do not affect the final converged ground state [Hj. 

For C2 = or c = 0, it is easy to check that SMA is always valid since the 
energy functional is symmetric with respect to spin component index. The fractional 
populations for each component is therefore the same as for a homogeneous system, i.e. 
given by ( ^-"°+'^ , uq, ^~"°~"^ ) if S = 0, and (m, 1 - m, 0) if S > 0. 

For ^^Rb and ^^Na condensates, which are believed to be ferromagnetic C2 < 
(c < 0) and anti-ferromagnetic C2 > (c > 0) respectively. Figure El gives typical 
density distributions of spacial mode function, p(r) = |0j(r)p. Both sub-figures clearly 
indicate that SMA is no longer valid. To get an overall idea of the validity of SMA 
we plot in Figure El the overlap integrals of our mode functions with respect to the 
SMA mode function 0sMA(r) as determined from a scalar GP equation with a nonlinear 
coefficient oc Cq (due to the symmetric Hs only) For a '^^Rb condensate, we see 

the overlap is close to unity when B is small, therefore, SMA remains approximately 
applicable. But it becomes increasingly bad with the increase of B. We thus conclude 
that the SMA remains reasonable in a weak magnetic field while it is clearly invalid in a 
strong 5-field. In fact, our numerical results confirm that the stronger the i?-field, the 
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Figure 7. The same as in Fig. El but now comparing the spin asymmetric energy 
Ea = C2(F)^/2 — (rjo + 'r]){Pz) + H^z) '^it^^ ^^'^ spin symmetric one Hs. 



worse the SMA gets. For typical system parameters, the dividing hne occurs at a S-field 
of a fraction of a Gauss when the system magnetization is not too small or too large. 
For a condensate with anti-ferromagnetic interactions, it was found earlier that SMA is 
violated in the limit of both large N and A4 even without an external 5-field, while the 
case of = presents an exception where SMA remains strictly valid for B = |14j . 
Figure ini shows the overlap integral for a ^^Na condensate, indeed we see SMA is invalid 
except at = where all atoms are in the |0) component. Remarkably, despite the 
seemingly large deviations from the SMA (as in Fig. the spin asymmetric energy 
term remains very small in comparison to the spin symmetric term as evidenced in Fig. 

13 

Figure |H1 shows the dependence of fractional populations on the fractional 
magnetization for a ^^Rb (left column) and a ^^Na condensate (right column) at 
different S-fields. For ^^Rb atoms, these curves resemble the same dependence as for 
a homogeneous system where SMA is strictly valid. Nevertheless, we find the densities 
of mode functions can become quite different, i.e. SMA is not valid in general. For 




Figure 8. Fractional population for each spin component of a *^Rb (left column) and 
a ^'^Na (right column) condensate. The values of B-field from top row to bottom are 
B — 0,0.1, 1.0 (Gauss). The atomic parameters are the same as in Fig. 5. The solid 
lines with plus signs denote the |+) component, the lines with triangles are for the 
|— ) component and the lines with open circles for the |0) component. The vertical 
dashed lines in (e) and (f) indicate the critical value rric, the boundary between the 
two distinct regions discussed in the text. In (d), rric = 0. 



^^Na atoms, the fractional component populations at different 5-fields again follow the 
results as obtained previously for the homogeneous case. When B = [as in Fig. 
ISfd)], the ground state distribution clearly obeys the same earlier (homogeneous) result 
(^^, 0, ^^), including the special case when = where it becomes (^"^, uq, ^"^) 
with ^0 G [0, 1]. For B [as in Fig. |Hl^e) and (f)], our numerical solutions reveal 
again two distinct regions; one for m < rric where all three components coexist, and 
another one for m > rric where only two components (|+) and |— )) coexist. We find 
that rUc increases with the 5-field, and is of course limited to rUc < 1. We conclude 
that despite the fact a harmonic trap induces spatially inhomogeneous distribution to 
condensate density, thus breaks the SMA in general, the overall ground state properties 
as measured by the fractional component distributions follow closely the results as 
obtained previously for the homogeneous case. Physically, we believe the above results 
can be understood as fractional populations relate to integrals of wave functions over 
all spaces, during which differences between wave functions can be averaged out. When 
only the |+) and |— ) components coexist, in fact, the two constraints on N and Ai 
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Figure 9. The B field dependence of the critical fractional magnetization rric as 
computed numerically for a ^^Na condensate in a harmonic trap. The smooth curve 
corresponds to the result of S = 2c[l — -y/l — ] (as from the homogeneous case) with 
an appropriately adjusted coefHcient c (or density). The atomic parameters are the 
same as in Fig. 5. 




always give the fractional population n± = {1 ± m)/2 if A^^o = 0. 

5. Conclusion 

We have revisited the question of the mean field ground state structures of spin-1 
condensate in the presence of a uniform magnetic field. For a homogeneous system, 
when c = 0, there exists in general only two nonzero components |+) and |0), except 
when B = where the ground state solution becomes indefinite; for ferromagnetic 
interactions when c < 0, the ground state in general has three nonzero components; 
when c > as for anti-ferromagnetic interactions, except for m = 0, there are two 
regions: one for 6 > 2c[l — ^1 — ] where three nonzero components coexist and one 

for 6 < 2c[l — y'l — ml ] where only two components coexist. Inside a harmonic trap, 
these results remain largely true, although the SMA becomes generally invalid. We find 
interestingly (see Fig. the B field (or the 6) dependence of the critical value that 
separates the two and three component condensate regions, remains almost identical as 
that given by the analytical formulae 6 = 2c[l — y'l — ] for a homogeneous system. 
In a sense, this also points to the validity of the use of a mean field description, as the 
number of atoms is really large (10^). 
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